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ABSTRACT 

We model the surface brightness profiles of a sample of 119 Abell Brightest 
Cluster Galaxies (BCG) , finding a generalised deVaucouleurs R^/" law, where n 
is a free parameter, to be appropriate. Departures from the R^^^ law are shown 
to be a real feature of galaxy profiles, not due to observational errors or coupling 
of n with the other model parameters. BCG typically have values of n greater 
than 4. The shape parameter n is shown to correlate with effective half-light 
radius, such that the larger BCG have larger values of n. This continues a trend 
noticed amongst ordinary elliptical galaxies and dwarf ellipticals, such that the 
fainter galaxies have smaller values of n. 



Subject headings: galaxies: fundamental parameters — galaxies: structure 
galaxies: elliptical and lenticular, cD 



1. Introduction 

The structure of the first-ranked members of Abell clusters, or brightest cluster 
galaxies (BCG), has long been of interest given that these galaxies can be used as standard 
candles for exploring the large-scale structure of the universe ( [Humason, Mayall, fc Sandage 
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1956; [Sandage 1972ii [1972E| ; |Gunn fc Uke"T975| ; Lauer & Postman 1994). |Hoessel (1980" 



for example, developed the L — a relationship, which relates BCG structure to their 
luminosities through fixed physical apertures, thus reducing the cosmic scatter in the 
aperture luminosities. Understanding the L — a relationship has been of particular interest 
recently, given its use by Lauer & Postman ( |1992| , 1994) to probe the linearity of the Hubble 
fiow and to detect large scale bulk fiows out to the 15,000 km/s scale. The structure of 
BCG may also refiect galactic cannibalism (|Ostriker fc Tremaine 1975| ; [Hausman fc Ostriker 



19781 ; Hoessel 1980; [Schneider, Gunn, fc Hoessel 19831) , where other cluster members may 



be captured through dynamical friction, changing both the luminosity and the shape of the 
dominant BCG. 

While to first order BCG appear to be ordinary, if highly luminous, elliptical galaxies, 
in many ways they are a special class of objects. [I'remaine fc Kichstone (1977)| showed, for 
example, that their luminosity function was not consistent with them being drawn simply 
as the brightest member of a standard [Schechter (1976)1 luminosity function. Pemler (1976)1 
showed that as a class, BCG differed structurally from other giant ellipticals of similar 
luminosity, having more extended limiting radii at a given total luminosity than other giant 
elliptical galaxies, ^chombert (1986)1 conducted an extensive survey of BCG brightness 



profiles, finding them to be shallower than those for more ordinary elliptical galaxies, such 
that the BCG were always more extended at a given surface brightness level than would 
be suggested by simply scaling by total luminosity. IHoessel, Uegerle, fc Schneider (1957)1 
also concurred with this, showing BCG to follow a shallower (and tighter) relationship 
between effective radii. Re, and effective surface brightness, Je, than do ordinary elliptical 
galaxies QOegerle fc Hoessel 19911 ; Graham 1996). This distinction in form was additionally 



interesting as ^chombert (1986")1 showed it to arise in the relatively brighter mid-portions of 
the BCG profile, regardless of whether or not the BCG had additional envelopes at faint 
surface brightness making them cD galaxies as well. A cD galaxy is a giant elliptical that 
has a separate extended low surface-brightness envelope, which is evident as an infiection 
in the brightness profile (Oemler 1976) — typically at /iy ~ 24 or greater ( IKormendy fc 



Djorgovski 195^0 . 



Schombert (1986)' emphasized that the classic r^^^ law, introduced by de Vaucouleurs 
( 11948 , 19531) to describe elliptical galaxies at all luminosities, was often a poor match to 



the BCG profiles, in many cases fitting the profiles only over a restricted range of surface 
brightness. Looking at the profiles presented by ISchombert (1986]| it becomes apparent 



that many of the BCG would be better fitted by power laws rather than r^/^ laws. This is 
also found to be the case in Ledlow fc Owen's (1995) work with galaxies in rich clusters. 
This appears to have little to do with whether or not a BCG is also a cD galaxy. An 
additional caveat is that because a constant power law will rise above an r^/^ law at large 
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radii, a cD envelope may be erroneously detected as a separate component in r^/^ plots, 
even though a single power law could describe the BCG completely. In short, classification 
of a BCG as cD or not is problematic; our present investigation, however, is not affected by 
this issue as this apparently has little to do in any case with the distinct over all structural 
properties of BCG from those of giant ellipticals. 

If some BCG are better described by power laws, and others by r^/^ laws, an inspection 
of the profiles presented by |Schombert (1986]| shows that BCG fall along a continuum 



between the two forms. We are thus motivated to advance the [Sersic (1968)| form, which 
includes both laws as a better description and generalization of BCG brightness profiles. 
In the Sersic form, /(r) oc exp(— (r/i^g)^/"), where n is a free parameter. As n — 1 from 
the n = 4 de Vaucouleurs form, the Sersic form becomes increasingly exponential, while as 
n —>■ oo, the Sersic form becomes a pure power law. 

Young & Currie (1994)' fitted the Sersic form to a sample of dwarf elliptical (dE) 
galaxies, showing a correlation between galaxy absolute magnitude and the shape parameter 
n, such that the fainter galaxies had the smaller values of n. The same pattern in profile 
shape had also been noted by |Davies et al. (1988)| . Working with a sample of Fornax 
low surface brightness (LSB) diffuse dwarf elliptical and spheroidal galaxies, they showed 
a trend between n and the logarithm of the scale-radius Re such that the scale radius 
increased as n increased. Caon, Capaccioli, and D'Onofrio (1993) showed that this trend 
extended into the domain of ordinary elliptical galaxies. It has also recently been shown 
to hold for the bulges of spiral galaxies by [Andredakis, Peletier, fc Balcells (1995)| and 
Courteau, de Jong, fc Broeils (199F) 



We are now in an excellent position to explore the systematics of the BCG structural 
properties as [Postman fc Lauer (1995)| present a complete and volume-limited sample of 



BCG out to 15,000 km/s. The data was obtained with large-area CCDs and reduced with 
photometric uniformity in mind, given the goal to use this data set to detect the subtle 
effect of bulk flows on the L — a relationship residuals (Lauer & Postman 1994). Here we 
show that the correlation between n and scale-radius (Caon, Capaccioli, and D'Onofrio 
1993) extends to BCG, where n is typically greater than 4, whereas for the dE it is seen 
to be less than 4. The existence of this relation across such a large range of galaxy sizes 
must be telling us something fundamental about the structure and formation processes of 
elliptical galaxies in general. 

In the next section we present some of the theory behind the i?^/" model. Section 
3 presents the BCG surface brightness data and the parameterised model flts to it. In 
Section 4 we discuss the results of the model fltting and compare this in the context of work 
on other galaxies in Section 5. We investigate how n relates to other measures of galaxy 
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structure, both local, Section 6, and global, Section 7. Our conclusions are given in Section 
8. 



2. The law 

The law gives the observed galaxy intensity, I, as a function of radius such that 
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where le is the intensity at the radius rg. The constant b is chosen so that re becomes the 
radius enclosing half of the total light from the galaxy. This generalised deVaucouleurs law 
was introduced by [Sersic (1968)| and further developed by [Ciotti (1991)| , the deVaucouleurs 
law has n=4 and b=7.67. Note that the constant b in this formula is a function of n. The 
generalised expression can be written as 



2.56„ 
/n(10) 



(r/re) 



1/n 



(2) 



with Te the scale radius, /Iq is the central surface brightness and 6„ the function of n given 
below. Again we select 6„ such that is the radius enclosing half of the total light for the 
^1/" model. For n=4, the well known deVaucouleurs R^^^ formula, = /io + 8.33(r/re)^/^ 
is recovered. The luminosity interior to the radius r is given by 



L(r) = Ierl2T^n{l - e) 



\2n 



7(2n, x) 



(3) 



for galaxy ellipticity e, and where 7(2n,x) is the incomplete gamma function with 

X = bn{r/rey^'^, defined by 

7(2n,x) = r e-H''''-^dt. (4) 



Thus, the value of 6„ is such that T{2n) = 7(272, where F is the gamma function. As 
given by [Capaccioli (1989") , this can be well approximated by 6„ = 1.9992?t, — 0.3271. 

The differences between an i?^/" model and an R^^^ model are best described 
graphically. In Figure |l|, we plot these differences for n ranging from 1 to 10. The models 
have been constructed to have the same total luminosity, in units of /g. The abscissa is in 
units of the half-light radius of the R^^^ model. As n climbs greater than 4, the curvature 
in the profile is steadily removed until it approximates a straight line. As n exceeds 10, the 
differences in the profiles become less marked. For these high values of n, the the portion of 
the model profile that is matched by observations is well represented by a power law. It is 
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seen that for the range in radius that can be matched by observations, values of n < 4 lead 
to a hump in the profile indicating that the E^/^ profile is too faint for the central parts and 
too bright to match the inner and outer portions of the profile. The situation is reversed 
for n > 4, such that the brighter galaxies have more light at larger radii than the R^/'^ law 
allows. 



3. Observational surface brightness profiles and fitted models 

The galaxies studied here are the BCG of the volume-limited sample of the 119 
Abell clusters known to be within 15,000 km/s, which were selected by [Lauer fc Postman 



(1994)1 to define a large-scale inertial reference frame. Full details of the galaxy selection. 



observations, and data reduction are given in [Postman fc Lauer (1995)] , but we summarize 



them briefly here. BCG candidates were selected from sky-survey images, and imaged in 
the Kron-Cousins R band under photometric conditions with large-area CCD cameras at 
KPNO and CTIO. The distance indicator used by [Lauer fc Postman (1994)] is based on the 



integrated luminosity of the BCG within the central 10h~^ Kpc metric aperture, thus final 
selection of a given galaxy as "brightest" among rival candidates is by metric rather than 
total luminosity. On this note, [Postman fc Lauer (1995)] emphasize that a number of their 



BCG identifications differ from those made earlier by [Hoessel (1980 

Postman & Lauer (1995)' measured profiles of the BCG using the multi-isophote fitting 
algorithm of P^auer (1986)| . Many of the BCG are part of multiple-galaxy systems; the 
multi-isophote algorithm solves for the brightness distributions of all overlapping galaxies 
simultaneously. Compact galaxies, stars, CCD defects, and so on, can also be excluded 
from the isophote fitting algorithm. Final accuracy of the profiles is limited by photometric 
calibration at small radii, and sky subtraction errors at large radii. [Postman fc Lauei] 



(1995)1 used repeat observations of the same galaxies to show that the basic random error in 
the BCG aperture magnitudes is only 0.014 mag; the surface brightness values for isophotes 
within the metric radius will be accurate to a similar level. Errors in the photometry due 
to sky subtraction errors will begin to dominate outside the metric radius, but should still 
be at the few percent level at the limiting isophote brightness adopted here of /ij? = 23.5. 



Postman fc Lauer (1995)[ did sky subtraction using a mode-estimator in the image corners. 



Some BCG were still contributing significant amounts of light at the edges of even the large 
fields of the CCDs used, in which cases additional observations were obtained offset from 
the primary galaxy image to measure the sky at larger angular distances from the BCG. 
Lastly, we note that the galaxy profiles are presented as observed. No K or extinction 
corrections have been applied (although [Postman fc Lauer (1995)[ did apply them to derive 
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BCG metric absolute luminosities). 

The model profiles have been fitted to the data outside of the central 3 arc seconds, 
due to the possible influence of core structure that is separate to the outer galaxy profile. 
We also have not used the profile data fainter than 23.5 magnitudes to be sure we are not 
affected by sky subtraction errors. In addition, this level of truncation in the profile ensures 
that our results are not a product of the extended haloes or envelopes that cD galaxies are 
known to have. For many of the BCG we had multiple images and could directly compare 
different profiles of the same galaxy for agreement. The model parameters obtained for 
galaxies with multiple images are found to be in agreement with each other within the 
errors. 

Equation § has been fitted to the semi-major axis surface brightness profiles of the 
BCG, using a simple error weighting scheme based on the S/N of each data point, via 
standard non-linear least-squares to solve for the three unknowns. We give here the 
superscript n to the value of /io and derived from the i?^/" formula, to prevent confusion 
with the values derived from the R^^^ formula. Ax^ ellipsoids were computed around the 
best fitting parameters, /io,Tg, and n. This was done by moving through a fine 3D grid 
of values and computing the value of at each point, rg was converted to log re and the 
ellipsoids then projected onto the relevant 2D plane. 

In addition, we fitted for the standard R^^'^ formula, where the value of n is fixed at 4 
and one solves for /xq and rg. We also explored the use of a power law to describe the light 
profiles. The method of least-squares was used to fit the profiles to 

jj{r") = A + B\og{r"), (5) 

where A is referred to as the intercept and is approximately the central surface brightness, 
being the value at r = 1". B is the slope of this power law. We divide B by -2.5 so that it 
reflects the slope of the log(Intensity)-log(radius) profile. Future references will refer to this 
modified value of B as being the slope. For the very flattened profiles, n > 15, the profiles 
lacked any significant curvature and it was appropriate to fit a power law of the above form 
to the data. 



4. BCG profiles 

The various profiles, R^^^, R^^^ and a simple power law, have been fitted to the BCG 
surface brightness profiles and are displayed in Appendix A. The residuals of the data 
about these best fitting models is shown clearly in Figure ^ for a handful of galaxies, as are 
the measurement uncertainties, based on S/N measurements, associated with the data. In 
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general, one finds that the R^^'^ law has too much curvature to match the data, resulting 
in a negative bowl shaped residual profile. The opposite is found for the fitting of a power 
law, where one generally finds the profile data has some level of curvature that cannot 
be accounted for with a simple power law, resulting in a positive hump in the residual 
profile. The i?^/" model with its free shape parameter can account for the differing levels of 
curvature and thus provides the best fit to the data, ironing out the large scale departures 
seen in the above two models. 

Shown in Figure are the best values of n plotted against logr" for the BCG. Also 
shown are the projected Ax^ = 9.21 ellipses about these points, corresponding to a 99% 
confidence region ( [Press et al. 198(j| ). It is pointed out that the profiles have additional 
wiggles in them that are not accommodated for by the i?^/" law and these exist at a level 
greater than the observational errors. As a result, the reduced values for each profiles 
optimal fit is larger than 1. Given that we have not under-estimated our errors, as shown 
by the good agreement of repeat measurements, this would imply that the model being 
fitted is inadequate. It is true that there are features/wiggles in the profiles that are not 
explained by the R^^"' law, nor are they explained by the R^^^ law. However, the bulk shape 
of the profile is described better by an i?^/" law than an i?^/^ law, as indicated in Table 1 
which shows the reduced values for the models fitted to the data. This is illustrated in 
Figure ^, which has the reduced values of each model's optimal fit plotted against the 
value of n from the i?^/" profile fit. One can see some general trends present. As one would 
expect, the R^^'^ profile fits are best when the i?^/" fits have n=4. The power law fits are 
better for larger values of n, which is to be expected as larger values of n mean a profile 
with less curvature. The i?^/" fits are better than those of the other two models and the 
quality of the fit appears to be independent of n. As an alternative measure of the errors, 
we set the reduced = 1 for the optimal fit and computed the la ellipses normalised to 
this level. These are also shown in Figure 0. 

Seeing effects are not responsible for the observed correlation of n with radius. We 
explored excluding different inner portions of the light profile, ^aglia et al. (19931 showed 



that the effects of seeing on the photometric properties of elliptical galaxies can extend as 
far as 5 seeing discs. Not using the inner 7.5", the same general trend between n and log re 
is still obtained. 

Not surprisingly, the i?^/" model is the best performer, as it can represent the R^^^ 
profile when n=4 and it can approximate a power law for large values of n, and fit for 
profiles of intermediate type. This can be attributed to having 1 more free parameter than 
the R^^^ or power law. What is important, is that such a variety in profiles are real, as 
indicated by the error ellipses in Figure ^. To explore this further, we simulated a pure 
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law profile out to IRg and added to it random gaussian noise, with standard deviations 
varying as a function of radius and being derived from the mean S/N errors of the 119 
BCG. We ran a Monte Carlo simulation, solving for r^,fiQ, and n each time with a new set 
of random noise to see if we could explain the range of values in r" and n observed with the 
real data. Figure || shows the entire cloud of solutions, not just the 1 or 2 ex elhpses. Clearly 
observational noise cannot explain the claimed trend of n with , suggesting that it is 
physically real and galaxies do indeed exhibit such a range in profile shapes. The possibility, 
mentioned by Kormendy (|1980| , |1982| ), that such a correlation may not be physically 
significant but rather due to coupling of the parameters in the fitting formula does not 
fully explain the observed variation in n, or the correlation of Ve with n. This analysis does 
however exclude the influence of sky subtraction errors, which are estimated to be at a level 
of 0.3%, and should not be significant. In addition, the S/N weighting scheme, employed in 
the model fitting, will down play possible contributions from sky subtraction errors. 

In fitting the i?^/" formula, it was not possible to obtain meaningful results for all 
profiles. The flattened profiles had large values of n, and implausibly large half-light radii. 
As the shape parameter increased, it lost its sensitivity and for values of n greater than 
about 15, the profiles resembled straight lines. That is, the profiles were better described 
by a power law. For about 40% of the data it was not appropriate to fit an i?^/" law, but 
rather a power law. This can be compared with the work of Ledlow fc Owen (1995)| who 
found for a control sample of 50 non-radio selected galaxies in rich clusters, 39% preferred 
a power law fit. The slope and intercept of the best fitting power law for each BCG is 
plotted in Figure ^ Not surprisingly, there is a slight trend for galaxies with a steeper 
profile to have a brighter central flux, represented by the /i intercept at r=l". As the power 
law is appropriate for BCG with n > 15 these galaxies have been marked with a star so as 
to distinguish them from the other galaxies that are not well represented by a power law. 
Both types occupy similar regions, with the BCG well fitted by a power law not occupying 
any special domain of this diagram. Of course the real galaxy profile must turn over and 
be truncated at some radius, otherwise it would be of infinite size and luminosity. But for 
the range in radius explored, some 40% of the BCG had surface brightness profiles of a 
power law nature. Ordinary elliptical galaxies and certainly dE do not show this feature 
over the same range in radius. It is pointed out that for large values of n, the half-light 
radius obtained from the i?^/" model is not really physical, but is however an expression of 
the slope of the galaxy profile measured by the data. 

We find no correlation between n and ellipticity. Following Uaon et al. (1993) , we have 
plotted n against the maximum ellipticity of the galaxy in Figure |^. We also used the 
ellipticity at the semi-major axis radius of 10 Kpc and again found no evidence for a 
correlation. |Caon et al. (1993]] found a range of ellipticities for n < 4 but for n > 4 the 
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ellipticity was typically less than 0.3. A possible reason for this different result could be 
that we used BCG whilst |Caon et al. (1993]| used E/SO galaxies, indicating a difference 



between the galaxy classes. Although this seems hard to understand, we note that |(Jaon et 



al. (1993") only had a dozen galaxies with n > 4 and only two with tt, > 10. They also found 



that this trend disappeared when they used the value of n constructed from the semi-minor 
axis profile. 

Figure || shows n plotted against the metric magnitude, being that enclosed by a 
circular aperture of 12.5 Kpc (we used a value of Hq = 80 km s~^ Mpc~^, and use the CMB 
reference frame here). A linear correlation coefficient of -0.17, at a significance level of 84%, 
suggest that there is little correlation between these two values for our BCG sample. This 
is not surprising given the small scatter in metric magnitudes for the BCG, being only 0.33 
magnitudes. That the BCG metric magnitude is not significantly infiuenced by the global 
galaxy profile shape, represented by n, eliminates the chances of another hidden variable 
giving scatter to this magnitude which has been used as a distance indicator by [Lauer fc 
Postman (1994)| and [Colless (1995)| in their studies of peculiar velocity fiows. 



5. i?^/" Universality 

What infiuences the profile shape of galaxies? Is it because of environment, changes 
arising due to dynamical evolution, differing formation history, or individual galaxy 
peculiarities such as dust clouds or rings, lenses, shells, ripples, etc.? The trend of n with 
galaxy size argues against galaxy peculiarities being responsible, as this would require such 
features to be correlated with galaxy size. The lack of a disk for the BCG eliminates the 
option that a disk is the responsible feature in determining the galaxy shape. However, 
the small wiggles in the profiles that are not accommodated by the i?^/" law are at a 
significance level greater than accounted for by the observational uncertainty and may be 
signs of such features. A close analysis of the two dimensional residual maps (images with 
the best fitting i?^/" law subtracted) could reveal such features. 

To investigate possible environmental effects on the BCG profile shape, we checked for 
a correlation between n and Richness Class, RC, ( |Abell 1958[ ) and n and Bautz-Morgan 



type, BM, ( |Bautz fc Morgan 19701 ). We only used those galaxies which had values of 
n < 15, being some 60% of our sample. No obvious trend was found for either case, as can 
be seen in Figure || and Figure respectively. We point out that the profile shape has been 
determined from the inner surface photometry (< 23.5 mag) and is thus free from outer 
galaxy distortions such as envelopes, which may be infiuenced by the environment. [Einasto 



|fc Caon (1993)1 have shown n not to correlate with galaxian density for their sample of 
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E/SO galaxies; therefore further restricting the possibihty of environment being responsible 
for the shape parameter n. 

BCG typically like values of n greater than 4. This is not only interesting in itself, but 
becomes more so when we note that studies of dE galaxies show them to have values of n 



typically less than 4. |Caon et al. (1993]| show that ordinary elliptical and SO galaxies show 
a range in n above and below 4. Continuing figure 5 of |Caon et al. (1993)| we add our data 
to produce Figure |Tl], showing the continued relationship between galaxy structure and size 
for dwarf galaxies, ordinary E/SO galaxies and BCG. |Caon et al. (1993)] fitted the i?^/" law 
to the B-band photometry of 33 E/SO galaxies. Whilst we used the Kron-Cousins R-band 
to image our sample of BCG, changes in profile shape due to color gradients are expected to 
be lost in the scatter of Figure |TT]. We note that the values of from |Caon et al. (1993}| are 



not those from their model but are derived directly from the galaxy light profile. The dwarf 



galaxies are a sample of 187 Fornax cluster dwarfs from Pavies et al. (1988)| . They fitted a 
variant form of the i?^/" law to the B-band images, such that /(r) = lQexp[—{r/A)^]. By 
setting N=l/n, one has RJ' = A"{2/N - 0.327)^/^. We used a value of Hq = 80 km s"^ 
Mpc~^ throughout in determining the value of Re Kpc. 

We thus have a parameter that traces structural differences amongst galaxies having a 
large range in size. The fact that the i?^/" profile is applicable over such a range (6 orders 
of magnitude in the mass) suggests that there are similar formation mechanisms present 
for all these galaxies. Computer modelling of the physics of violent relaxation has shown 
that the shape parameter n is dependent on the central potential of the models ( [Hjorth S 



Madsen 1995| ). Galaxies with larger potentials result in galaxies with larger values of n, 



and those with smaller potentials have values of n less than 4, in agreement with what the 
observations reveal. This suggests that the bulk shape to an elliptical galaxies luminosity 
profile is not due to its individual genetic peculiarities or its environment but rather 
something intrinsic to each galaxy, being its mass. This has previously been suggested by 
Young fc Currie (1994)| and [Andredakis, Peletier, fc Balcells (1995)| . We hope to investigate 



this further when we have velocity dispersion data for the BCG. This information can be 
used to estimate the mass of each galaxy ([van Albada, Bertin, fc Stiavelli 1995| , [Michard 



1980|), which we can then check for a correlation with n. 



6. Local profile structure a 

The shape parameter, n, describes the overall global shape of the galaxy light profile. 
There is another measure of a galaxies light profile curvature, referred to as the structure 
parameter by Poessel (1980)| and given by a = dlnL/dlnr\j.^. It is a measure of the galaxy 
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profile at a given sampling radius, r^. Lauer fc Postman (1994) and CoUess (1995) have used 



a because of its correlation with metric magnitude and hence ability to create improved 
standard candles for distance measurements. In each case, a was measured directly from 
the light profiles at a sampling radius of 10h~^ Kpc. It is, however, possible for one to 
compute a from the fitted model profile [(Graham 1996)| . Whether this is preferred or not. 



one gains insight into how a varies for different values of n. Working from equation ^ and 
assuming zero ellipticity for simplicity, it can be shown that 

a = — 77, (6) 

wy^Zn, x) 

A set of i?^/" profiles were constructed and a computed as a function of radius for each. 



The results of which are shown in Figure |12| For a fixed sampling radius, r^, a is seen to 



increase as the effective half-light radius increases. Working against this is the fact that as 
galaxies get bigger, ie. larger i?e, they tend to be described by a profile with a larger value 
of n, as seen in Figure ^ Now profiles with a larger values of n, are seen to have smaller 
values of a for the same ratio of sampling radius to effective half-light radius. The degree 
to which a changes for a given change in sampling radius or change of n is dependent on 
the part of the profile one is looking at. 

We have plotted the values of n against a for our sample of BCG that have values of 
n < 15 in Figure |l^. Values of a have been taken from [Lauer fc Postman (1994)| rather than 



computed from the fitted profile model. It seems likely that the reason no trend is evident 
between these two parameters is because of the competing situation described above. 



7. Another shape parameter 

An alternative approach to quantifying the systematic deviations from an R^^^ law 
was taken by Purkert (1993)| . The observational data in the range 0.6 < x < 1.1 was fitted 
by the form = /io + bx, where x = {r/X^y^^ and is the effective half-light radius 
of the model. In general, /io and b are dependent on the fitted range in x, which is in turn 
dependent on the value of X^. As we don't have a prior knowledge of the value of Xg, 
this dependence means that our expression must be solved via an iterative process until 
convergence, such that X"'^"' = X°''^(8. 3268/6)^, with b the best fitting slope using Xf'^. 



Burkert observed that many galaxy profiles had a characteristic dip/hump about the 
best fitting R^^'^ law. When the dip {fXgaiaxy — R^^^) was a minimum, it was found to turn 
over near x=0.8, and when it was a maximum it would be near x=0.9. See [Burkert (1993) 



for a fuller description and profiles. The above mentioned range in x was then cut into two 



- 12 - 



parts, with the division at Xcut = 0.8 if the dip was a minimum, or at Xcut = 0.9 if it was a 
maximum. 

Cdu Ou \ 

^ Ug[x.„,;1.i] Ue[o.6;.w] 1 /8.3268 (7) 

was shown to correlate with both absolute magnitude and the mean deviation of the galaxy 
data about the best fitting R^^^ law. Burkert showed that the brighter galaxies had the 
more negative values of 6b, and the fainter galaxies had the more positive values. 

This parameter reflects the curvature present in the galaxy profiles that the R^^'^ law 
can't accommodate. However, it's extraction is somewhat fiddly and requires the fitting 
function be applied several times, whereas the i?^/" profile is fitted once. A series of i?^/" 
profiles were created, and the method of Burkert applied to obtain the parameter 6b for 
each. The relation between n and 6b is shown in Figure |14[ The sense and amplitude of n 
and 6b are in agreement, such that the bigger brighter galaxies with n > 4 relate to 56 < 
and the galaxies with n < 4 have 6b > 0. 



8. Conclusions 

The nature of the i?^/" law is such that as n increases the galaxy surface brightness 
profile flattens, and for ri > 15 it is very well represented by a power law. BCG typically 
have values of n greater than 4 - that is their light profiles are less curved than the classic 
i?^/^ law; 40% of our sample is well described by a power law. The range in profile shapes is 
real and not due to noise in the galaxy profiles or due to a couphng of the three parameters 
in the model. There is a trend between n and the half-light radius, such that the larger 
galaxies have larger values of n. This trend appears to be a continuation of that noticed 
for dE galaxies, LSB galaxies, through to normal E and SO galaxies and on to BCG, 
suggesting some common physical processes must be at play in the formation of all of these 
galaxies. This global shape parameter, n, is shown to be independent of Richness Class and 
Bautz-Morgan type, suggesting that the galaxy environment (in so far as RC and BM type 
represent this) is not responsible for the shaping of the bulk distribution of stars in the 
galaxy. We note that our analysis excludes the outer envelopes of the cD galaxies. 

Whilst n is a global measure of the galaxies light profile, a is a measure the galaxies 
structure at a point, a is shown to be related to n and R^ in an opposing manner. As 
one moves to larger galaxies, n increases causing a to decrease but at the same time Re 
increases causing a to increase. The dominating factor depends on which part of the profile 
one samples the value of a at. 
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A. Brightest Cluster Galeixy Surface Brightness Profiles 

The surface brightness profile data for all 119 BCG, together with the best fitting E}/^ 
(dotted line), i?^/" (sohd line), and power law (dashed line) are displayed. The models were 
fitted to the error weighted data brighter than 23.5 magnitudes (Kron-Cousins R band) and 
outside 3 arcseconds (filled circles). Data outside this range being given by an open circle. 
The S/N errors are shown, but usually are contained within the data point symbol. The 
surface brightness level at 22.0 magnitudes is indicated for each profile by a small arrow on 
the right hand margin. No i?^/" profile is shown when n > 15. 

REFERENCES 

Abell, G. O. 1958, ApJS, 3, 211 

Andredakis, Y. C, Peletier, R. F., Balcells, M. 1995, MNRAS, 275, 874 
Bautz, L. P., & Morgan, W. W. 1970, ApJ, 162, L149 

Binney, J. 1989, in The World of Galaxies, ed. H.G.Corwin and L.Bottinelli, Springer- Verlag, 
Berlin, p.208 

Burkert, A. 1993, A&A, 278, 23 

Capaccioli, M. 1989, in The World of Galaxies, ed. H. G. Corwin, L. Bottinelli (Berlin: 
Springer- Verlag) , 208 

Capaccioh, M., Caon, N., & D'Onofrio M. 1993, MNRAS, 259, 323 

Caon, N., CapaccioU, M., & D'Onofrio, M. 1993, MNRAS, 265, 1013 

Ciotti, L. 1991, A&A, 249, 99 

CoUess, M. 1995, AJ 109, 1937 

Courteau, S., de Jong, R. S., & Broeils, A. H. 1996, ApJL, 457, 1 

Davies, J. I., Phillips, S., Cawson, M. G. M., Disney M. J., & Kibblewhite, E. J. 1988, 
MNRAS, 232, 239 



-14- 



de Vaucouleurs, G. 1948, Ann d'Ap., 11, 247 

de Vaucouleurs, G. 1953, MNRAS, 113, 134 

de Zeeuw, T., & Pranx, M. 1991, Ann Rev A&A, 29, 239 

Einasto, M., k Caon, N. 1993, MNRAS, submitted 

Graham, A. W. 1996, ApJ, in press 

Gunn, J. E., & Oke, J. B. 1975, ApJ, 195, 255 

Humason, M. L., Mayall, N. V., & Sandage, A. R. 1956, AJ, 61, 97 

Hausman, M. A., & Ostriker, J. R 1978, ApJ, 224, 320 

Hjorth, J., & Madsen, J. 1995, ApJ, 445, 55 

Hoessel, J. G. 1980, ApJ, 241, 493 

Hoessel, J. G., Oegerle, W. R., & Schneider, D. P. 1987, AJ, 94, 1111 

Kormendy, J. 1980, in ESO Workshop on Two-Dimensional Photometry, ed. P. Crane, K. 
Kjar (Garching:ESO), 191 

Kormendy, J. 1982, in Morphology and Dynamics of Galaxies, ed. L. Martinet, M. Mayor 
(Sauverny: Geneva Obs.), 113 

Kormendy, J., & Djorgovski, S. 1989, ARA&A, 27, 235 

Lauer, T. R. 1986, ApJ, 311, 34 

Lauer, T. R., & Postman, M. 1992, ApJ, 400, L47 

Lauer, T. R., & Postman, M. 1994, ApJ, 425, 418 

Ledlow, M. J., & Owen, F. N. 1995, AJ, 110, 1959 

Michard, R. 1980, A&A, 91, 122 

Oegerle, W. R., & Hoessel, J. G. 1991, ApJ, 375, 15 

Oemler, A. 1976, ApJ, 209, 693 

Ostriker, J. P., & Tremaine, S. D. 1975, ApJ, 202, L113 
Postman, M., & Lauer, T. R. 1995, ApJ, 440, 28 



- 15 - 



Press, W. H., Flannery, B. P., Teukolsky, S. A., & Vetterling, W. T. 1986, Numerical 
Recipes (Cambridge: Cambridge University Press) 

Sandage, A. 1972a, ApJ, 173, 485 

Sandage, A. 1972b, ApJ, 178, 1 

Saglia, R. P., Bertschinger, E., Baggley, C, Burstein, D., CoUess, M., Davies, R. L., 
McMahan, R. K., & Wegner, G. 1993, MNRAS, 264, 961 

Schechter, P. 1976, ApJ, 203, 297 

Schneider, D. P., Gunn, J. E., & Hoessel, J. G. 1983, ApJ, 268, 476 
Schombert, J. M. 1986, ApJS, 60, 603 

Sersic, J.-L. 1968, Atlas de Galaxias Australes (Cordoba: Observatorio Astronomico) 

Tremaine, S. D., & Richstone, D. O. 1977, ApJ, 212, 311 

van Albada, T. S., Bertin, G., & Stiavelli, M. 1995, MNRAS, 276, 1255 

Young, C. K., & Currie, M. J. 1994, MNRAS, 268, Lll 



This preprint was prepared with the AAS IM^jX macros v4.0. 



-16- 



Fig. 1. — a) Surface brightness profiles of an i?^/" law, for integer values of n ranging from 1 
to 10. b) The difference between the surface brightness profiles of an i?^/" and an 7?^/^ law, 
for n ranging from 1 to 10 in steps of 1. All profiles give the same total magnitude, in units 
of 4" 

Fig. 2. — We show the surface brightness profiles for a sample of BCG, together with the 
best fitting R^^^. (dotted line) i?^/" (solid line) and power law (dashed line). The residuals 
of the profiles about these models are displayed in a magnified portion below the profiles 
themselves, being the circles, stars and triangles for the respective models. 

Fig. 3. — upper panel) Plot of n versus logr^ (Kpc) for our sample of galaxies. (40% of the 
BCG had values of n > 15). middle panel) 99% confidence regions, lower panel) la ellipses 
after the reduced has been normalised to 1. 

Fig. 4. — Reduced values from the optimal fit of the i?^/"^, power law, and i?^/" models 
to the BCG profiles are shown here against the value of n from the optimal i?^/" fit. 

Fig. 5. — Cloud of solutions from a Monte Carlo investigation of the contribution to the 
n — log Tg relation from observational errors. 

Fig. 6. — Plot of the best fitting power law slope and intercept at r = 1" for the BCG. 
Those galaxies fitted with a value of n > 15 from the R^^^ model have been marked with a 
star, those fitted with a value of n < 15 are marked with a square. The la error bars being 
smaller than the graph markers. 

Fig. 7. — Plot of n against the maximum ellipticity of the BCG. 

Fig. 8. — Plot of n versus the metric magnitude contained with 10h~^ Kpc. 

Fig. 9. — The galaxy shape parameter, n, is plotted against the galaxy cluster Richness 
Class, RC. Also shown are the mean values for each class and the associated la confidence 
interval. 

Fig. 10. — The galaxy shape parameter, n, is plotted against the Bautz-Morgan, BM, type 
for the galaxy cluster. Also shown are the mean values for each BM type and the associated 
la confidence interval. 

Fig. 11. — The galaxy shape parameter, n, is plotted against the logarithm of the effective 
half-light radius. Re Kpc, for a sample of dwarf galaxies (Davies et al. 1988), triangles, 
ordinary E/SO galaxies (Caon et al. 1993), filled circles, and our BCG data, open circles. 

Fig. 12. — The structure parameter, a, is plotted as a function of radius for different values 
of the shape parameter, n, from the i?^/" model, r^/rg is the ratio of the sampling radius, 
at which a is computed, to the effective half-light radius of each respective model. 

Fig. 13. — Plot of n versus a, where a is taken from Lauer & Postman (1994). 
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1.02 


0.53 


20.95 


1.17 


17.33 


-2.02 


2.01 


3581 










1.25 


22.96 


3.34 


16.84 


-1.43 


2.42 


3656 










1.17 


22.12 


4.34 


15.97 


-1.47 


2.57 


3676 


1.6 


0.86 


20.88 


3.24 


1.06 


21.91 


3.58 


16.92 


-1.60 


3.85 


3677 


7.3 


1.11 


22.95 


1.18 


0.92 


22.01 


1.55 


17.56 


-1.74 


1.64 


3698 


11.2 


1.29 


23.47 


2.56 


0.79 


21.10 


3.41 


15.88 


-1.69 


3.00 


3716 


10.1 


2.73 


28.51 


2.19 


1.62 


23.94 


2.64 


17.87 


-1.30 


2.44 


3733 










1.30 


23.20 


3.35 


17.58 


-1.46 


2.46 


3736 


7.4 


1.98 


25.05 


1.67 


1.51 


23.00 


2.68 


17.06 


-1.44 


2.81 


3742 


12.2 


1.05 


22.74 


2.86 


0.65 


20.77 


3.76 


15.56 


-1.78 


3.28 


3744 


12.8 


1.38 


23.82 


2.13 


0.88 


21.46 


2.68 


16.76 


-1.75 


2.31 


3747 


5.4 


1.26 


22.88 


2.44 


1.11 


22.18 


2.86 


16.74 


-1.54 


3.57 


3869 


6.2 


0.98 


22.04 


1.51 


0.88 


21.53 


2.03 


16.93 


-1.78 


2.45 


4038 


10.9 


1.79 


25.50 


2.19 


1.08 


22.26 


3.17 


16.78 


-1.57 


2.77 


4049 










1.13 


22.41 


4.16 


16.72 


-1.58 


3.25 


4059 


4.7 


1.92 


24.52 


2.25 


1.78 


23.93 


2.37 


17.63 


-1.24 


3.29 



NOTE. — The reduced value of the best fitting R^^'^ , R^^^ , and power law model is shown for each BCG, along with the parameters 
from the fit. Re and R^ are the effective half-light radii of the R^^^ law and the R}^'^ law respectively, with jLtg and jLt^ the surface 
brightness at these points. For values of n > 15, the R^^'^ fitting function fails to produce meaningful results. 
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